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We propose a new type of Josephson vortex ratchet. In this system a Josephson vortex moves 
in a periodic asymmetric potential under the action of a deterministic or random force with zero 
time average. For some implementations the amplitude of the potential can be controlled during 
the experiment, thus, allowing to tune the performance of the system and build rocking as well as 
flashing ratchets. We present a model describing the dynamics of the fluxon in such a system, show 
numerical simulation results, and discuss the differences between conventional and Josephson vortex 
ratchets. The investigation of this system may lead to the development of the fluxon rectifier — a 
device which produces a quantized dc voltage from colored noise (non-equilibrium fluctuations). 
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I. INTRODUCTION 



To extract useful work from random motion was a dream of mankind since the days when the Brownian motion was 
recognized. Unfortunately, the second law of thermodynamics forbids to extract energy "for free" from equilibrium 
thermal fluctuations (white noise), which was didactically demonstrated by R. Feynman in his Lectures Never- 
theless, one can extract useful work from non- equilibrium or time- correlated (colored) noise "not paying" for it, using 
so-called ratchets, i.e., systems with an asymmetric periodic potential Recently there was a boost of activity in 
this field related to the experimental investigation of directed motion in biological systems, so-called Brownian motors 
which, e.g., move muscles or transport vesicles in a cell ||. In the latter case the probable mechanism of operation is 
the motion of kinesin molecules along the surface of microtubules, which can be mapped to the motion of Brownian 
particles along a one-dimensional ratchet potential with the period 8.2 nm Q. The non-equilibrium energy is supplied 
by chemical reaction of splitting of adenosine triphosphate (ATP) which takes place close to the kinesin molecule. 

In addition to the application of ratchets as noise rectifiers, it was suggested to use them for very efficient separation 
of small objects with different mobility, e.g., DNA molecules, viruses, etc. H§. The particle separation is based on 
so-called deterministic ratchets , where the particles move in a certain direction under the action of a deterministic 
force with zero time average. Moreover, changing the force profile one can reverse the direction of the particle motion 
jjj . The classification and discussion of different types of ratchet systems can be found in Ref . . 

In this paper we focus on Josephson ratchets which are of particular interest because (a) the directed motion results 
in a dc voltage according to the Josephson relation and (b) these systems can operate at very high frequencies up 
to about 100 GHz. As a first example we mention the asymmetric dc SQUID where the equation of motion for the 
Josephson phase (difference of quantum mechanical phases) corresponds to the motion of an imaginary particle in a 
2D ratchet potential. Such SQUID ratchets have been proposed (|] and studied experimentally Another type 
of Josephson ratchet investigated recently is a ID array of Josephson junctions with spatially modulated properties 
(□J. A Josephson kink (vortex), which can move along the array, can be considered as quasi-particle in a ID ratchet 
potential. 

Here we propose a new class of Josephson ratchets which further develops the idea of a kink in a ID array. The 
proposed system consists of a ID long Josephson junction (LJJ) which may be bent in the a6-plane [see Fig. [l](a)] or 
have variable width w(x). Here and below x is a curvelinear coordinate along the junction. The fluxon (Josephson 
vortex) moving along the junction, from a mathematical point of view, is a topological soliton. It has its own mass, 
velocity, and other particle-like properties |L2| . We study the motion of a fluxon in L J Js in a ratchet potential which 
can be formed either by applying an external magnetic field and bending the junction properly or by modulating its 
width w(x). To provide the required periodicity of the potential, the junction is topologically closed in a loop. Such 
a geometry is similar to the well known annular Josephson junction E|-[l7j in which a fluxon moves in a sinusoidal 
potential created by a magnetic field. Using a more elaborated shape JTa ] one can form an asymmetric potential with 
the possibility to control its amplitude by changing the amplitude of the external magnetic field. An alternative idea 
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of magnetic field modulation using a specially shaped control line is proposed in Ref . |f{| . The directional motion 
of a fluxon can be detected by measuring the dc voltage across the junction which is, due to the Josephson relation, 
proportional to the average velocity of a fluxon. 

Before discussing fluxon dynamics in a ratchet potential we would like to stress the difference between conventional 
ratchets and Josephson vortex ratchets. First, the fluxon, although it has a topological charge and in general behaves 
like a particle, is a nonlinear wave, i.e., it can change its shape rather strongly as well as emit electromagnetic waves. 
Second, the fluxon dynamics is usually studied in the underdamped limit which is opposite to the overdamped case 
which was considered for the majority of work on Brownian particles in a ratchet potential. Small damping may result 
in chaotic dynamics even in the deterministic case and, in the case of the fluxon, even without any potential 
Thus, to have a well defined behavior of a fluxon ratchet, one has to work in the overdamped or in the weakly 
underdamped limit. If one uses conventional Nb-AlO^-Nb technology to fabricate the LJJ, this requirement means 
that the working temperature should be very close to T c . As an alternative one can use junctions with intrinsically 
high damping such as SINIS LJJs Q or high-T c LJJs technology Q which allows to fabricate LJJs of required 
topology. Third, if we consider multi-particle dynamics, the strong repelling interaction between fluxons plays an 
important role and must be taken into account. 

This paper is organized as follows. In section [n] we derive the equations for the dynamics of the Josephson phase 
in a bent LJJ of variable width w(x) in the external magnetic field. We also discuss different kinds of fluxon ratchets, 



their advantages and drawbacks. The numerical simulation results are presented in III. Section IV concludes this 
work. 



II. THE MODEL 



Here we derive the generalized perturbed sine-Gordon equation which takes into account the curvature of the LJJ 
in the afo-plane [see Fig. 0(a)], the uniform magnetic field H applied in the plane of the junction (in b direction) as 
well as the modulation of the LJJ width w(x) along its length x. 

We start from a discrete representation of the LJJ shown schematically in Fig. |f|(b) . The Kirchhoff equations for 
the Josephson phases in the cell and for the currents in one of the nodes are: 



(x + dx) — 4>(x) 
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where 4>(x) is the Josephson phase at point x of the junction, $(x) and & e (x) are the total and the external magnetic 
flux applied to the cell, respectively, L{x) is the inductance of the piece of the junction electrodes between x and 
x + dx, Il{x) is the current in the electrodes, i.e. through the inductance L(x), I e (x) is the externally applied bias 
current, and I(x,t) is the current through the Josephson junction. The particular expression for I(x,t) depends on 
the adopted JJ model and is introduced later. 

Assuming that the interval dx is infinitesimal, we can r ewrite Eqs. ([[]) and (g) in a differential form using the 
following expressions: 
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where Hod' is the inductance of one square of the superconducting electrodes 0], d! « 2Al is the effective magnetic 
thickness of the junction |24| , n is the unit vector normal to the plane of the junction cell as shown in Fig. [l](b) , A sa 2 Al 
is the effective penetration depth of the magnetic field into the junction |24| ], and Al is the London penetration depth 
of the superconducting electrode. We assume that the films are spatially uniform so that d! and A are independent 
on x. 

Substituting Eqs. into Eqs. (|l|) and (||) we can rewrite the latter in a differential form as 



9^ 
dx 



2tt 

d! L (x) 
dx 



HqH(x)A 



Hod' 



w(x) 

= w(x) [j e (x) -j(x)] 
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Excluding Il(x) from Eqs. (Q) and (|J), we get the equation which describes the dynamics of the Josephson phase 
in the system: 



$0 
Z7T 



(9) 



Here and below, the subscripts t and x, if any, denote the derivatives with respect to time t and coordinate x, 
respectively. Note, that we didn't include any particular model of JJ into our equation up to now, which is a definite 
advantage of this derivation procedure. In the case of the simple RSJ model, one should substitute j(x) which is the 
sum of the supercurrent, normal (quasiparticle) current and displacement current densities: 

j(x) = j c sin(0) + ^-<t> t + Cp-fa (10) 

ATTN, Z7T 

into Eq. (^). Here j c , R and C are the critical current density, specific resistance and specific capacitance of the 
junction, accordingly. In this case Eq. (^) can be rewritten in a form which resembles the usual sine-Gordon equation 
& 

\ 2 j<t> xx - Lof&t - sin(0) = uj% - l(x) + QH x {x) + [QH{x) - X 2 (/) x ] , (11) 

F w(x) 



where Aj = ^/(f>o/(27r/Ltoj c d') is the Josephson penetration depth, uo p — yj2irj c / ($oC) is the Josephson plasma 
frequency, tu c = 27rj c i?/$o is the so-called critical frequency, j(x) = j e (%)/jc is a normalized bias current density, and 
Q = 27r/x AA^/$o- 

For theoretical investigation of the system we introduce standard normalized units, i.e., we normalize the coordinate 
to the Josephson penetration depth Aj, and the time to the inverse plasma frequency w" 1 . After such simplifications, 
Eq. ( |TT| ) can be rewritten as |2^,|2(|: 

4>xx - 4>tt - sin(0) = a<j> t - j(x) + h x (x) H ^— [h(x) - (j> x ) , (12) 

w(x) 

with the damping coefficient a — uj p /uj c = l/^/~fT c , and the field h normalized in a usual way as 

tl c \ 

H c \ — $o/( 7r A t oAAj) is the first critical field (penetration field) for a LJJ which is, in fact, equal to the field in the 
center of the fluxon. The normalized velocity is given in natural units of cq — \jto p , where cq is the so-called Swihart 
velocity. The normalized voltage V = <pt is given in units of $o w p/(2tt). From now on all quantities are given in 
normalized units. 

In comparison with the usual perturbed sine-Gordon equation, Eq. ( |l2| ) contains 3 additional terms. The term 
h x (x) describes the effect of the applied magnetic field when the junction is bent in the afe-planc. The second term 
[w x (x)/w(x)] 4> x comes from the width modulation. The last term [w x (x)/w(x)] h(x) describes the mixture of both 
and does appear only when both, field modulation due to curvature and width modulation are present. 

It can be checked by a direct substitution into the Euler-Lagrange equation 

d_ dL_ + d_ dC_ _dC _ Q 

dt dcf> t dx d(j) x dej) 

that the Lagrangian density 

- [ ^" %)]2 -(l-cos^)l, (15) 




results in the equation of motion (|12j) without the acpt and 7 terms which describe dissipation and external force and 
therefore are not included in the Lagrangian density. Of course, the Lagrangian density (|l5|) can be obtained directly 
from Fig. |l|(b) and the RSJ model. From Eq. (|l|) one can see that instead of the usual potential energy term cf> x /2 we 
now have w{x) [4> x — h(x)] 2 /2, i.e., actually 3 terms. The first, w(x)(j) x /2, is the obvious generalization of the usual 
potential energy term to the case of variable width w(x). The second term w(x)h(x) 2 /2 is a constant term due to the 
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applied field and is not related to the fluxon motion or other Josephson phase activity in the junction. In fact, there 
are no traces of this term in Eq. The third term h(x)<p x w(x) represents the part of the potential energy density 

which we are interested in and which we are going to exploit to build a system with a ratchet potential. 
One of the solutions of the sine-Gordon Eq. (|12|) with zero r.h.s. is a soliton (fluxon) 

4>{x,xq) = 4 arctanexp(a; — xq), (16) 

with the center situated at point xq. We consider non-relativistic motion i.e., dxo/dt <§C 1. Further assuming that the 
fluxon profile (|l6| ) does not change much due to the r.h.s. of Eq. (|l2]) which acts as a perturbation, we can get the 
explicit expression for the potential energy U(xq) as a function of fluxon coordinate Xq. For this purpose, we note 
that in expression ( |l5| ) for the Lagrangian density the second and the third terms correspond to the potential energy 
density U(x,xo) (with opposite sign). The potential energy U(xq) is obtained by substituting <j>(x,xo) from Eq. (|lq ) 
into U(x,xq) and integrating over x. Thus, we get 



Aw{x) 2w(x)h(x) 



U(x ) = / 4^ S K \ dx. (17) 

cosh (z-xq) cosh(x-xo) 



' — oo 



The first term corresponds to the potential energy due to width modulation, the second due to shape, field and width. 

The first possibility to form a ratchet potential is to apply no magnetic field (h = 0) and to vary the width w(x) 
of the junction. In this case the potential energy will be given by the first term of Eq. (|l7|). Moreover, when the 
junction width w(x) does not change much over the distance comparable with the fluxon size, the potential 

U{x Q ) » 8w(x ) (18) 

just repeats the w(xq) profile. From Eq. (|l^) it follows that if the width changes as a saw-tooth, so does the potential 
energy, except for the vicinity of the saw-tooth's infinite slope where one has to calculate the convolution according 
to Eq. (0). Here and below when mentioning saw-tooth profile we mean a saw-tooth with finite positive slope and 
infinite negative slope. In the case of the saw-tooth where the width changes from wq to wo + Aw at the point xq = 
(infinite slope), the potential energy in the vicinity of the point xq = is U(x$) = 8wo + 4Au> [1 + tanh(xo)]. The 
practical implementation of such a ratchet would look like an annular LJJ with the outer edge having the shape of a 
circle and inner edge having the shape of one turn of a spiral. This geometry has the advantage that the corresponding 
potential can be made ideally saw-tooth- like except for the smearing due to the convolution in Eq. (|lj), which is a 
common feature of all fluxon based systems. In some sense, this system is an analog of the Josephson ratchets based 
on ID arrays []Tl"l| . Unfortunately both types of ratchets do not allow to control the potential height during experiment 
which can be considered as a disadvantage. 

The second possibility is to keep w constant, to apply a magnetic field and to bend the junction in the afo-plane. 
In this case the first term of Eq. ( |l7| ) gives a constant and we have to consider only the second term. Again, if h{x) 
changes slowly in comparison with the fluxon size, Eq. ( [l7| ) is simplified to 

U(x ) = -2irwh(x ). (19) 

In the well known case of the ring shaped junction, the field h(x) = hg cos(#) = ho cos(x/R) [8 is the angle between H 
and n, as shown in Fig. |lj , and therefore we get a symmetric potential. If we deform the ring properly, the potential 
can be made asymmetric as desired. Possible experimental shapes are shown in Fig. ^. The advantage of this kind 
of ratchets is that one can control the amplitude of the potential during experiment by varying the amplitude of 
the magnetic field. The possibility to tune the potential height allows to implement so-called flashing ratchets || by 
applying either an ac magnetic field using a coil or just placing an rf antenna close to the junction so that the H of 
the emitted electromagnetic wave will have non-zero component in the a6-plane. 

Note, that the terms j(x) and h x (x) in Eq. ( [l2| ) from a mathematical point of view play the same role. Therefore 
in the experiment the field h(x) can be substituted by a properly chosen additional bias current 7 P (a;) = — h x (x), 
which has zero average in space. The inverse mapping is also valid, but the bias current with non-zero average maps 
to a non-periodic field (potential) with linearly growing component, which does not belong to the class of ratchets. 



In section [II .A we, in fact, use (7) 7^ (these brackets denote spacial averaging), but only to test the asymmetry of 



the potential. When we demonstrate the real operation of the ratchet in the section [II. B, we have (7) = 



III. SIMULATION RESULTS 



In this section we study the fluxon dynamics in a Josephson ratchet of the second type (h ^ 0, w = const) for a 
saw-tooth field profile h(x). This ideal asymmetric profile is not only of academic interest because it should show 
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the best figures of merit for fluxon ratchets but also can be quite closely reproduced in a real system [see shape in 
Fig. 0(b) and corresponding h(x) in Fig. [j|(a)]. All simulations were performed using an explicit numerical scheme for 
Eq. ( |l2| ) using a LJJ of the normalized length i = 20, with damping coefficient a — 0.2 (weakly underdamped limit). 
The numerical technique and simulation software are discussed in detail elsewhere fjPijl . 



A. Probing the asymmetry of the potential 

First, we probe the asymmetry of the potential by calculating critical current 7 C vs. potential height fto for the 
case of one trapped fluxon [7^[(fto)] and for the case of no trapped fluxons [7^0(^0)]- The superscripts '+' or '— ' 
correspond to opposite directions of applied bias current. 

On the basis of the model derived in the previous section we can understand how these dependences should look 
like for arbitrary magnetic field profile h(x). In general, due to the loop- like geometry, the left and right tails of a 
fluxon can interact. We assume that the junction is long enough and this interaction in negligible. This situation 
is equivalent to the long periodic system where the fluxons are separated by a large distance I 3> Xj. Furthermore, 
we represent the field h(x) as hoTC(x) and map it to the equivalent additional bias current 7 P (x) = —hoH x (x). An 
example of h(x) and 7 P (x) derived from the geometry in Fig. ||(b) is shown in Fig. [|(a). 

If there is a fluxon in the junction, 7(2;) and 7n(a;) translate into a driving force ^(a^o) and a potential (pinning) 
force F p (xo), respectively, acting on the fluxon |2^] 

F 7 (x ) = f +0 ° 2 J. (X)W , dx = 2ir 1W , (20) 
J_ OQ cosh(a;-xo) 



F p (x )= ( + °° y )W dx = h a f(x )w. (21) 
J^oo cosh(x - XQ) 

The corresponding potential U(xq) and force f(xo) are shown in Fig. |^(b). For the sake of simplicity we suppose that 
7 does not depend on x, but our discussion can be easily generalized to the case when the r.h.s. of Eq. (po| ) is equal to 
2t:T(xo)w. The pinning force ( pl| ) can be also obtained directly from expression ( |l7| ) for the potential energy. When 
we increase 7 the fluxon is pinned while these two forces can compensate each other, i.e., Fj(xq) + F p (xq) = 0. The 
depinning happens for 

7 = 7ci = -fto/(*i)/27r, (22) 
where Xi is the coordinate at which f(xi) has a minimum [see Fig. ||(b)]. We assume that 7 > and f(xi) < 0. Thus 



from Eq. (22) we see that 7 c i(fto) looks like a straight line which starts from the origin. 

On the other hand, regardless of the presence of a fluxon, the Josephson phase evolves under the action of the 
net current 7 + 7 P . From Fig. ^(a) it is clear that at point x 2 , where "t P (x) is maximum, this sum can exceed 1 
for some value of 7 = "f c and the junction switches to the resistive state. Using the 7 p = —h$H x (x) the condition 
7c + -i p {x 2 ) = l gives 

7c(fto) = 1 + h H x (x 2 ) = l-ho \H x {x 2 )\ , (23) 

i.e., also a straight line, but with the negative slope. 

Summarizing, j c i(ho) has two branches (|22| ) and (p3|) , and follows the one with the lowest critical current for given 
ho- For fields < ho < ft.*, j c i{ho) is given by Eq. Q22|), and for ho > ft* by Eq. (p3[), where the field 

K = f(x 1 )/2n + n x (x2Y (24) 

is the field where these two dependences intersect. The dependence 7 c o(fto) has only one branch (^) since the fluxon 
depinning mechanism is absent. 

When one applies the current 7 in the opposite direction, the dependences qualitatively look the same but the 
particular values of slopes can be different in the case of an asymmetric potential, if 7^(2^2) 7^ —7^(0:4) and f(x\) ^ 
—f(xa) (see Fig. ||). Thus, the measurement of 7^(fto) and 7^(fto) gives direct information about the asymmetry of 
the field and of the potential, respectively. Note that an asymmetric potential may have 7^(2:2) = —7^(2:4) and/or 
f{x\) — —f(xa), i.e., it does not reveal its asymmetry in the measurements of 7^o(fto) and 7 c i(fto)- Inversely, the 
asymmetry of 7^(fto) and 7^(fto) is a clear indication of asymmetry in the system. 
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We stress that the first branch of 7^(/io) an d the second branch of 7^(/io) show the asymmetry of the field h(x), 
while the first branch of 7^(/io) shows the asymmetry of potential U(xq), i.e., already after convolution with <f> x . 
While h x (xo) may have appreciable asymmetry, U(xq) may be almost symmetric. 

We should note that the above analysis is valid for an arbitrary field profile h(x) as long as h x (x) can be considered 
as a perturbation in the r.h.s. of Eq. (|l2|). However, for the case of the saw-tooth potential this is not the case, 
since h x {x) has a <5(x)-like behavior at the point of the junction where the saw-tooth has infinite slope. Therefore our 
analysis is valid only qualitatively. In particular the dependences j^(ho) still have linear slopes but the values of the 
slopes can not be calculated using simple integrals as shown above because the fluxon shape will differ from ( fl6|) quite 
considerably, as we saw in numerical simulations. Therefore we present the results of direct numerical simulation of 
Eq.@. 

The simulated 7^(ft-o) an( f 7^i(^o) f° r & n ideal saw-tooth field profile h(x) are presented in Fig. H. The magnetic 
field is given in terms of the force amplitude h§2-K j i. If the system is closed in a loop with a normalized circumference 
I, the h(x) can be written in the form h^TL{2irx/ 1) where Ti. is geometry dependent 2n periodic function normalized 
to 1. Therefore the term h x (x) in Eq. ( fL2| ) is 

h x {x) = ho^H' (j-x\ , (25) 

where H'(£) = dH(£)/d£ is a periodic function. This means that the potential force scales inversely proportional to 
the length of the junction I. So, to get rid of the length dependence and make our results valid for any length I 3> 1, 
we present our results as a function of ho2n/£. 

From Fig. ^ one can see that the depinning current 7ci(0) = and then grows linearly with field, as expected from 
the theory, up to some field h*2Tr/£ w 1.8. After that it decreases linearly with field, also according to our prediction, 
and, finally, exhibits multiple branches corresponding to the multi-fluxon states. 

As it should be in the case of a saw-tooth potential, the slopes of 7^(/io) and 7^(/iq) are not symmetric, and the 
ratio of the slopes is about 4. Ideally, if the fluxon were a real particle we could expect an infinite force necessary 
to push the fluxon out of the well in the direction against the infinite slope and, therefore, an infinite ratio of slopes 
on jfi(ho) curves. In practice, due to the convolution (|l~7|), we have a finite force and a finite ratio of slopes on the 
7^(/io) dependences. The finite ratio of slopes on jf Q (h ) as well as of the second slopes 7^(/io > ^*) is observed 
because perturbation theory does not apply as mentioned above. 

The next step is the understanding of the dynamics of the fluxon in an asymmetric potential. One of the simplest 
observations which can be made numerically as well as experimentally is the examination of the fluxon trapping current 
7t r , i.e., the minimum current at which a fluxon still moves along the system, not being trapped by the potential. 
Obviously, in the underdamped system the trapping of the fluxon by the potential will not take place while the kinetic 
energy of the fluxon exceeds the height of the potential, and one should not see any difference in the fluxon trapping 
currents ^(ho) for opposite bias current directions. In the strongly overdamped case the fluxon dissipates energy so 
fast that to move further under the action of the driving force, the driving force should always overcome the maximum 
value of the potential (trapping) force. In an asymmetric potential the maximum force created by the potential is 
different for opposite directions of fluxon motion and we should expect a difference in the fluxon trapping current. 
The simulated ~f^.(ho) curves for a saw-tooth potential are shown in Fig. ^. Even in the slightly underdamped case 
a = 0.2 the dependences for positive and negative direction of fluxon motion (bias current) differ quite considerably. 
We also note that 7^(/io) almost coincides with 7^(^-0) (see Fig. ^) for ho < 2 which means that for positive direction 
of bias there is no hysteresis on the I-V characteristic when the system switches from zero voltage state to the state 
with moving fluxon and back. Instead, for the negative bias the hysteresis is present as can be seen in Fig. || 

B. Motion due to monochromatic force 

In the following we investigate the motion of a fluxon in a system with the saw-tooth field profile h{x) and w = const 
under the action of a monochromatic ac bias current 7 = 7 ac sin(ojt) in Eq. (p"2"|). If as a result of an ac drive the 
fluxon starts to move around the junction, we can estimate the characteristic frequency of this process. It is clear that 
the maximum velocity of the fluxon is equal to the Swihart velocity (1 in normalized units) therefore the maximum 
revolution frequency is w ma x = 2n/£ w 0.314. 

In the limit of low frequency (quasi-static) drive u> -C w max we can calculate the net average velocity of the fluxon, 
just integrating the current-velocity characteristics ^(7) (see Fig. ||) at a given value of field ho as: 

1 [ 2 * 

u(7ac) = tt I w(7 ac sm(r)) dr, (26) 
2tt Jo 
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where the integration path should follow the proper hysteretic branches of the u(j) dependence. 

In the case when lo is comparable with the maximum frequency of the fluxon rotation around the junction w max , 
we perform direct numerical simulations. The dependences u(j ac ) calculated numerically for different values of u> are 
shown in Fig. 0(a)-(c) for ho2n/£ — 0.5, 1.0 and 2.0, respectively. The upper curves marked as "qs" (quasi-static) are 
calculated using Eq. (ptf). All other curves for ui — 0.01-0.1 were obtained by means of direct numerical simulation 
of Eq. (Il2|). Such values of frequencies were chosen taking into account that they should be less but comparable with 
ii) ma . As we saw in simulations, at u> > 0.2 [in the case of Fig. |?](a) already for ui — 0.1] the rectification is suppressed 
almost completely. 

From Fig. [?] one can see that the dependence u(7 ac ) for quasi-static case is piecewise with jumps at the values of 
bias current where switching between different hysteretic branches occurs. For example, let us focus on ho2w/£ = 1.0 
see Fig. 0(b)]. In this case the switch between different hysteretic branches occurs at 7 which is equal to 0.522, 0.696, 
0.916 (c.f. Figs. ||, ^ and [?]). Several characteristic regions can be distinguished. First, in the region of small values of 
7 ac , where u — 0, the amplitude of the ac current is not big enough to exceed the pinning current and the fluxon is 
localized in the well. At higher ac bias up to 7 ac = 0.522 the positive dcpinning current is exceeded and the system 
switches to the fluxon step, i.e., the fluxon escapes from the well and starts to move around the junction. Note that 
during a period of ac drive the fluxon can only make integer number of turns and at the end of the period will be 
again localized in the well. This results in the quantization of the average velocity u, as can be seen as steps on all 
curves except for the quasi-static one. The step number corresponds to the number of turns done by the fluxon during 
one period of ac drive. The voltage rectifier based on this principle will give a quantized voltage V n — n& u; /2tt 
(in physical units) which is defined by fundamental constants and the applied frequency. The corresponding velocity 
quanta in Fig. are given by u n = nuj£/{2ir). This gives ui « 0.0318 for £ = 20 and to = 0.01. Note that the voltage 
quantum increases with u> and for u> = 0.01 five quanta (turns per ac cycle) give the same voltage as one quantum for 
uj = 0.05. For the quasi-static case the voltage quantum is infinitesimal so that the quasi-static curves look smooth. 

A further increase of the ac bias amplitude 0.522 < 7 ac < 0.696 results in a decrease of the dc voltage because 
during the negative semi-period of the ac drive the fluxon starts moving in the opposite direction so that the efficiency 
of the fluxon ratchet drops. 

The big peak at 0.696 < 7 ac < 0.916 in Fig. ^(b) is related to the switching of the junction to the resistive state 
during the positive semi-period while during the negative semi-period the junction stays in the zero-voltage state or 
on the fluxon step. When the amplitude of ac drive gets big enough to switch the system to the resistive state also 
during negative semi-periods, the dc voltage drops again and in the limit of very strong ac drive approaches zero. A 
large amplitude of ac bias 7 ac 3> 1 implies that the system spends almost the whole period in the positive or negative 
resistive state and only a tiny fraction of the period on the asymmetric part of the current-velocity characteristic at 
low currents, so that the resulting average velocity is close to zero. In addition, such a strong ac drive results in 
chaotic dynamics, as can be seen in Fig. [7| for large values of 7 ac . 

Figure also shows that the rectification effect decreases when the driving frequency increases and approaches the 
maximum fluxon rotation frequency w max . 

We point out that the performance of the fluxon rectifier depends on the chosen potential amplitude ho- For 
comparison in Fig. ^(a)-(c) we show the rectification characteristics for different values of magnetic field amplitude 
2-Khty/£ = 0.5, 1.0 and 2.0 (see Fig. [|). We see that the potential depth affects the amplitude of ac drive at which 
rectification appears, at which it(7 ac ) has maximum and it affects the relative location of the two regions corresponding 
to the fluxon rectification regime and resistive rectification regime. We note, that the resistive rectification regime in 
principle gives a larger effect which is inversely proportional to the damping coefficient a, i.e., it can be made rather 
big. However, this is a rather trivial effect which can be obtained using any non-symmetric I— V curve. The fluxon 
rectification regime, which is the main subject of this study, gives smaller but a quantized voltages. This can be an 
advantage for some applications. 

IV. CONCLUSION 

We proposed a new type of Josephson vortex ratchets, where the motion of a fluxon along a long Josephson junction 
closed in a loop can be considered as the motion of a quasi-particle in a ratchet potential. The derived model suggests 
several implementations of the fluxon ratchets, and, in particular, the one where the amplitude of the potential 
can be controlled during experiment. Since the fluxon is a soliton (non- linear wave) moving in the underdamped 
medium with asymmetric potential, we expect some non-trivial effects related to (a) wave properties of the fluxon, 
(b) underdamped dynamics and (c) interaction between fluxons. As a first step, we performed a numerical analysis 
of the effective potential seen by a fluxon, and checked the rectification of a monochromatic signal. We found that 
a fluxon rectifier produces a quantized voltage, with a quantum given (due to the Josephson relation) only by the 
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fundamental constant $o and the driving frequency to. The experimental investigation of the proposed system based 
on (Nb-Al-A10a;)-Nb Josephson junction technology is in progress. 
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FIG. 1. Piece of LJJ (di and d,2 are the thicknesses of the superconducting electrodes and di is the thickness of the 
insulating tunnel barrier) which is bent in the afc-plane and has variable width w(x): (a) 3D view of geometry, (b) schematic 
representation using discrete elements. 

FIG. 2. Possible shapes of LJJs (top view) that provide a ratchet potential when the magnetic field H is applied in the 
direction shown by the arrow. 

FIG. 3. Illustration to the explanation of the 7^)(feo) and 7^(/io) dependences, (a) h(x) and 7p(x). (b) U(xq) and f(xo). 
These curves correspond to the LJJ shape shown in Fig. ^(b). 

FIG. 4. Normalized critical currents J^o(ho) (no trapped fluxons) and 7^(fro) (one trapped fluxon) vs. magnetic field 
amplitude ho for the saw-tooth magnetic field profile h(x) in the annular LJJ of length I — 20. Dotted lines show the values of 
magnetic field at which rectification shown in Fig. R(a)-(c) was calculated. 



FIG. 5. Normalized fluxon trapping current 7t r (^o) vs. magnetic field amplitude ho for the saw-tooth magnetic field h(x) 
in the annular LJJ of t — 20 and a = 0.2. Dotted lines show the values of magnetic field at which rectification shown in 
Fig. (?](a)-(c) was calculated. 

FIG. 6. Current-voltage characteristic of LJJ with one trapped fluxon shown as the dependence of fluxon velocity 
u = Vdc^/27r on the amplitude of the bias current 7 obtained for a — 0.2. 

FIG. 7. The dependence of the dc voltage across the junction Vdc given in terms of average fluxon velocity u = Vdc^/27r on 
the amplitude of the ac drive 7 ac for different normalized frequencies ui — qs, 0.01, . . . ,0.1 of the ac drive and a = 0.2. The 
curves corresponding to the different frequencies are intentionally shifted by 0.0636 (two velocity quanta at u) = 0.01) relative 
to each other for the sake of better visibility. The insets of (a) and (b) show magnified views of the fluxon rectification region. 
Figs, (a), (b) and (c) were obtained for 2nho/l = 0.5, 1.0 and 2.0, respectively. 
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